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Relying on the quantization rule of Raab and Friedrich [Phys. Rev. A (2009) in press], we derive 
simple and accurate formulae for the number of rotational states supported by a weakly bound 
vibrational level of a diatomic molecular ion. We also provide analytic estimates of the rotational 
constants of any such levels up to threshold for dissociation and obtain a criterion for determining 
whether a given weakly bound vibrational level is rotationless. The results depend solely on the 
long-range part of the molecular potential. 
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I. INTRODUCTION 

Molecular systems bound by a potential which varies 
asymptotically with the inverse power of the distance, r, 
between their constituents 



V(r) 



c 



with n > 2 



(i) 



are amenable to an accurate analytic semiclassical 
(WKB) treatment, as long as the system's states are suffi- 
ciently ensconced within the potential energy well. How- 
ever, for states near threshold for dissociation, the WKB 
approximation fails, as the system's classical action, pro- 
portional to momentum, no longer exceeds Planck's con- 
stant. And yet, it is the near-threshold states that have 
come to the fore recently, through the work in cold-atom 
physics where such states arise in photo- and magneto- 
association [1, 2, 3] or other types of "assisted" collisions 
of (ultra)cold atoms [4] or atomic ions [5, 6, 7, 8]. There- 
fore, a considerable effort has been expended at amending 
the WKB approximation to also allow for tackling near- 
threshold states analytically. A leading approach is that 
of H. Fricdrich et ai, who showed that, firstly, the bound- 
state eiegenencrgies, Eb, can be expressed in terms of the 
quantization function, F(Eb) [9], which relates the state's 
integral quantum number v to the generally non-intergral 
quantum number, vth, of a state exactly at threshold, via 



F(E b 



vth - v 



(2) 



Secondly, they were able to find an explicit analytic form 
of the quantization function for attractive inverse-power 
potentials with n = 6 [10, 11], and, most recently, with 
n = 4 [12]. Note that the binding energy, Eb = D — E v , 
with D the dissociation energy and E v the energy of the 
vibrational level v, is thus positive, Eb > 0. 

In our previous work [13], we have shown that for each 
vibrational level, v, the rotational angular momentum, 
J, can take a critical value, J*, such that the vibrational 
level is pushed up to threshold, thereby causing the level's 
binding energy to vanish. Hence the angular momentum 
J in excess of J*, J > J*, dissociates the molecule, cf. 
Figure 1. Furthermore, we have shown that the critical 



angular momentum is related to the quantization func- 
tion by 



J* = F(E b )(n - 2) 



(3) 



The corollary of Eq. (3) is that the number of rotational 
states supported by a weakly bound vibrational level of a 
dimer is given by the integer part of the critical angular 
momentum, J max = Int[J*]. By making use of the ex- 
plicit form of the quantization function of refs. [10, 11] 
for n = 6, we were able to evaluate J* and estimate the 
rotational constant B for near-threshold states of 85 Rb2. 

Here we take advantage of the recently derived explicit 
form of the quantization function for the n = 4 long- 
range potential and analyze the rotational structure of 
highly-excited H^~ and 133 Cs^~ molecular ions. The accu- 
racy of the approach based on the quantization function 
is compared with that of an exact numerical solution of 
the corresponding Schrddinger equation. 



II. ROTATIONAL STATES OF WEAKLY 
BOUND MOLECULAR IONS 

The quantization function for the n = 4 case takes the 
form: 

F{E b ) = F th {n) + F W ( K ) [>„.(«) + F WKB (k)] (4) 

where the individual terms, herein introduced for conve- 
nience, comprise a near-threshold dependence, 

2bn - (pn) 2 
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an "interpolation" term, which gives a smooth transition 
between 1ow-k and high-K behavior, 



F ip ( K ) = 



(G 6 nf + (G 7 k) 7 
1 + (G 6 k) 6 + (G 7 k) 7 ' 
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a term which corrects the reflection phase due to the 
potential of Eq. (1), 



F cr {n) 
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TABLE I: The w, 
Eqs. (5) and (6) 
G 7 = 1.338059. 



coefficients of Eq. (7). The coefficients of 
ire: b = 1, p 2 = G 6 = 1.622576 and 
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and a pure WKB contribution, 



WKB 



(«) 



(8) 



cf. Eq. (38) of rcf. [12]. In Eqs. (4)-(8), the dimensionless 
wavenumbcr k is defined by 



k = k 



/C 4 2m 
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(9) 



with k = y/2mEb/hthe wave vector and m the diatomic's 
reduced mass. The parameters b, p and u in Eq. (4) are 
listed in Table I (cf. Tables I, II of rcf. [12]). Note that 
in order to avoid confusion with the rotational constant 
(defined below), we changed the symbols Bq^, used in 
ref. [12] for the adjustable-length parameters, to Gq^. 

The quantization function (4) can be simplified by ne- 
glecting the F t h(n) and F cr (n) terms and setting Fi p (n) 
to ~ 1, which results in the Leroy-Bernstein (LB) ap- 
proximation [14]. Neglecting only the F t h{n) and Fi p (n) 
terms yields the improved Leroy-Berstcin approximation 
(iLB), which accounts for short-range deviations of the 
true potential from V(r) of Eq. (1), see refs. [10, 15]. The 
assumptions about the various terms of Eq. (4) inherent 
to the approximations are listed in Table II. 

Neglecting any coupling of the molecular rotation, we 
can estimate the rotational constant, B, from the ro- 
tational energy, BJ*(J* + 1), required to promote the 
vibrational level bound by E to threshold 



B 



E h 



J*(J* + 1) 



(10) 



The values of the rotational constant B obtained from 
Eq. (10) for 133 CsJ and H+ are listed in Tables III and 
IV together with the essentially exact values, B exact . The 
latter were calculated from 



TABLE II: Terms of the quantization function of Raab and 
Friedrich (RF), Eq. (4), inherent to the LB and the iLB 
approximations. Also shown are the ranges of the reduced 
wavenumber n wherein the approximations apply. See text. 
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be derived from Eq. (3) with J* 
for a level to be rotationless is 

E b < d 4 



mV2 C 4 1/4 



1 . Hence the criterion 



(12) 



The parameter has the same value for all potential 
wells with an 1/r 4 tail, namely d± — 2.9105. Within 
the LB and iLB approximations, the coeeficients d 4 LB = 



2.9096 and d 



LB 



0.7386. These values come close to 



those given in Table III of ref. [13]. 



III. RESULTS AND DISCUSSION 
A. 133 Cs ^ molecular ion 

As an example, we analyzed the rotational structure 
of the 133 Cs 2 ion, for which an accurate potential energy 
curve was published recently [8]. Table III compares, for 
the last three vibrational levels of 133 CsJ, the values of 
the critical rotational angular momentum J* and rota- 
tional constant B, calculated from Eqs. (3) and (10), with 
exact results. The tabic also lists the values obtained by 
the Lcroy-Bcrnstcin and improved LeRoy-Bcrnstein ap- 
proximations. One can see that for all the states consid- 
ered, the predicted values of J* come close to the exact 
values, as do the values of J* LB - However, J£ B , given by 
the purely semiclassical LB theory, are quite off the exact 
values. One can also see from Table III that the model's 
estimate of the rotational constants B is within 25% of 
the exact value. In the case of the LB approximation, the 
deviation of J£ B from the exact value happens to be in 
the direction such that it improves the agreement of the 
LB rotational constant with the exact one; however, this 
improved agreement has to be regarded as serendipitous. 
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with the vibrational wavefunctions obtained from a nu- 
merical solution of the Schrodinger equation for the po- 
tentials of refs. [8] and [16] . 

The value of the binding energy Ef, for which a vibra- 
tional state is unable to support molecular rotation can 



H, molecular ion 



We have also looked at the other end of the mass scale 
and tackled the rotational structure of the H^ ion. The 
exact values of J* and B were obtained by numerically 
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FIG. 1: (Color online) A schematic illustrating the role of the centrifugal term in the effective potential, V e g(r) = — + 

^7^73-^- The energy splittings have been exaggerated for clarity. Shown is the position of a rotationless vibrational level, 
v(J = 0) (dashed line), as well as the position of the same level when pushed up by the centrifugal term to threshold, 
v(J = J*) = vth(J = 0) (full line, at threshold). When the rotational angular momentum J exceeds the critical value J*, the 
centrifugal term pushes the vibrational level above threshold, v(J > J*) (full line, above threshold), thus leading to dissociation. 



TABLE III: Critical angular momenta J* and rotational con- 
stants B obtained for the last three vibrational states of the 
133 Cs 2 molecular ion in different approximations; Eb and B 
are given in 10 -8 cm -1 . See also Table II and text. 



V Eb J* Jexact J*LB JIb B B exa ct BiLB BlB 

371 5.7 0.63 0.63 0.63 1.01 5.53 4.02 5.51 2.81 

370 499.4 2.63 2.63 2.63 3.08 52.4 35.5 52.4 39.7 

369 3731.1 4.62 4.64 4.62 5.10 143.5 96.4 143.5 120.0 



TABLE IV: Critical angular momenta J* and rotational con- 
stants B obtained for the three least-bound states of the Hj" 
molecular ion in different approximations; Eb and B are given 
in cm -1 . See also Table II and text. 



V E b J* Jtxact J*LB JIb B B exact BiLB BlB 

19 0.707 1.27 1.33 1.27 1.70 0.24 0.15 0.24 0.15 

18 23.411 3.60 3.70 3.60 4.07 1.41 0.92 1.41 1.14 

17 153.67 6.03 6.56 6.03 6.51 3.62 2.10 3.62 3.14 



solving the the Schrodingcr equation for the potential of 
rcf. [16]. Table IV compares the exact critical angular 
momenta and rotational constants with those obtained 
from the models. Although the potential well of the hy- 
drogen molecular ion is rather shallow, the near-threshold 
rotational structure is still governed by the long-range 
potential tail, as attested by a reasonable agreement of 
both J* and B with the exact results. Since none of the 
last vibrational states lies in a "pure near-threshold re- 



gion," the LB and iLB approximations are not too far off 
cither. 

IV. CONCLUSIONS 

The advent of the physics of translationally cold atoms 
and molecules has secured a new prominence for long 
range interactions [17]. The r -4 attraction features most 
notably in the ion-induced dipolc interaction between 
an ion and and an atom or molecule, as well as in the 
Casimir-Polder potential of an atom interacting with a 
surface [18]. The bound states near threshold are of 
particular interest. This is because such states are non- 
classical, often spending most of their lifetime in the far- 
out, classically forbidden region of the potential [4]. Also, 
molecular species, whether ionic or neutral, in such states 
can be probed using nonrcsonant laser light, by "shak- 
ing" [19]. On the other hand, threshold states deter- 
mine low-energy scattering behavior, subject to studies 
in traps [5]. 

Within this context (and beyond), it is helpful to pos- 
sess the means the assess the rotational structure near 
threshold. Herein, we provide simple formulae that en- 
able such an assessment to be carried out. These for- 
mulae yield the critical value of the rotational angular 
momentum needed for dissociation; the rotational con- 
stant; and a criterion for determining whether a given 
vibrational state can support molecular rotation. 
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